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Dissociation  energies  for  the  step-wise  removal  of  hydrogen  from  methanol,  CH3OH  ->  CH3O  ->  CH^O  —  CHO  -» CO, 
are  obtained  by  many-body  perturbation  theory.  The  heat  of  formation  of  CH3O  is  predicted  as  2  ±  3  kcal  mole  '.  Com¬ 
puted  geometries  are  in  excellent  agreement  with  experiment,  where  available,  and  provide  a  prediction  for  CH3O. 
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1.  Introduction 

A  major  task  for  theoretical  chemists  is  the  devel¬ 
opment  of  methods  to  predict  energy  differences 
with  chemical  accuracy.  Most  quantum  chemists 
agree  that  accurate  prediction  of  relative  energies 
requires  application  of  theories  that  include  electron 
correlation  effects,  effects  not  treated  in  self-con¬ 
sistent-field  (SCF)  calculations  [1-4] .  Estimates  of 
molecular  correlation  energy  have  been  obtained  by 
using  the  method  of  configuration  interaction  to 
improve  the  wavefunction  obtained  from  an  ab  initio 
SCF  calculation.  Recently,  however,  theoretical 
techniques  based  upon  the  linked-diagram  theorem 
have  been  extended  to  molecular  applications  [1,2, 
4],  and  these  methods  provide  an  alternative  ap¬ 
proach  to  the  correlation  problem.  Several  compari¬ 
sons  of  the  various  theoretical  methods  have  ap¬ 
peared  [1,2,4]. 

This  report  describes  the  application  of  many-body 
perturbation  theory  (MBPT)  to  predict  the  dissocia¬ 
tion  energies  for  a  series  of  reactions  that  occur  in 
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flame  zone  combustion  models.  The  reaction  set 
comprises  the  unimolecular  reactions  that  oxidize 
methanol  to  carbon  monoxide: 


CH3OH  -*■  CH3O  +  H,  (1) 

CH30  -*  CH20  +  H,  (2) 

CH20-*CH0+H,  (3) 

CHO  -*  CO  +  H.  (4) 


Reactions  (1)  and  (2)  involve  the  methoxy  radical. 
Since  heat  of  formation  data  are  available  for  the 
other  moieties  in  those  two  reactions,  accurate  pre¬ 
dictions  of  the  dissociation  energies  for  these  reac¬ 
tions  will  enable  us  to  predict  the  heat  of  formation 
for  the  methoxy  radical. 

The  theoretical  techniques  used  are  outlined  in 
section  2,  and  the  results  for  the  series  of  dissociation 
reactions  follow  in  section  3. 


2.  Outline  of  the  computations 

For  the  closed-shell  molecules  included  in  this 
study,  a  restricted  Hartree-Fock  (RHF)  wavefunction 
represents  the  ground  state,  whereas  unrestricted 
Hartree-Fock  (UHF)  wavefunctions  represent  the 
ground  state  for  open-shell  molecules.  The  integrals 
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over  atomic  functions  were  computed  using 
Dunning’s  4$3p  contraction  of  Huzinaga’s  9s5p 
gaussian  primitive  set  for  first-row  atoms,  and 
Dunnings  scaled  (5  =  1.2)  3s  contraction  of 
Huzinaga’s  4s  primitive  set  for  hydrogen  (5,6) .  A 
single  set  of  d-type  polarization  functions  (gaussian) 
augments  the  atomic  basis  sets  for  carbon  and  oxygen, 
with  exponents  0.75  and  0.85  respectively  [7] .  A  set 
of  p-type  polarization  functions  (gaussian),  a  =  1.0, 
augments  the  hydrogen  basis  set.  All  integral  calcula¬ 
tions  were  performed  using  the  MOLECULE  integral 
program  [8] .  Electronic  structure  calculations  were 
performed  using  the  GRNFNC  and  UMBPT  pro¬ 
grams  \ 

The  total  energy  reported  for  each  molecule  per¬ 
tains  to  structural  parameters  that  have  been  opti- 

*  The  program  GRNFNC,  written  by  G.D.  Purvis,  does  SCF 
iterations  and  integral  transformations.  The  program 
UMBPT,  written  by  R.J.  Bartlett  and  G.D.  Purvis,  does 
MBPT,  CCD  and  VP  Cl. 


Fig.  1.  Definition  of  bond  lengths  and  of  bond  angles  for 
CH3OH,  CII3O,  ClljO,  and  CllO. 


mizcd  relative  to  double-excitation,  frozen-core, 
fourth-order  MBPT  calculations.  The  structural 
parameters  optimized  for  each  molecule  are  defined 
in  fig.  1.  Computed  structural  parameters  are  com¬ 
pared  with  available  experimental  values  in  table  1. 

In  the  absence  of  experiment,  we  offer  a  prediction 
of  the  methoxy  structure  which  we  estimate  to  be 
accurate  to  0.01  A  in  bond  length  and  2—3°  in  bond 
angles. 

We  report  dissociation  energy  predictions  for 
three  levels  of  calculation:  SCF,  fourth-order  MBPT 
including  all  double-excitation  diagrams,  D  MBPT(4), 
and  fourth-order  MBPT  including  all  single-,  double-, 
and  quadruple-excitation  diagrams  that  contribute 
at  that  order,  SDQ  MBPT(4). 

The  prediction  of  dissociation  energies  requires 


Table  1 

Optimized  structural  parameters  for  molecules  included  in 
this  study  a) 


Parameter 

This  work 

Experiment 

CH3OHb> 

n 

1.093 

1.0937 

'1 

1.093 

1.0937 

a 

107°  9' 

108°  32' 

methyl  tilt 

2°  12' 

30  ,2. 

1.432 

1.4214 

'1 

0.963 

0.963 

P 

ior  2' 

108°  2' 

CH3O 

'1 

1.085 

r 2 

1.081 

108°  42' 

<*2 

112°  12' 

Pi 

103°  48' 

Pi 

1 1 1°  30' 

S| 

1.405 

CHjOc> 

T\ 

1.102 

1.099 

J] 

1.211 

1.203 

a 

116°  11’ 

116°  30' 

CHOd) 

r  1 

1.111 

1.125 

1.188 

1.175 

a 

124° 

124°  57' 

COe> 

*i 

2.125 

2.127 

a)  Bond  lengths  in  A. 

Experimental  values:  ref.  (12). 

c)  Experimental  values:  ref.  (13). 

d)  Experimental  values:  ref.  [14]. 

e)  Experimental  values:  ref.  [15]. 
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accurate  computation  of  energy  differences.  In  addi¬ 
tion  to  the  change  in  the  electronic  energy,  vibra¬ 
tional,  rotational,  and  translational  energy  differences 
must  also  be  obtained.  We  assume  that  rotational  and 
translational  energy  effects  are  well  approximated  by 
the  equipartition  of  energy  among  those  degrees  of 
freedom.  Since  vibrational  frequencies  are  established 
for  all  the  molecular  species  addressed  in  this  study 
except  the  methoxy  radical,  empirical  zero-point 
energy  changes  are  used  for  these  cases. 

Characterization  of  the  methoxy  radical  was  the 
most  difficult  problem  encountered  in  this  study. 

Bent  et  al.  [9]  obtained  optimal  structural  parameters 
for  the  Jahn— Teller  distorted  ground-state  conforma¬ 
tion  after  extensive  D  MBPT(4)  calculations.  An  anal¬ 
ysis  of  the  potential  energy  function  obtained  in  that 
study  provided  guidelines  to  estimating  the  vibra¬ 
tional  frequencies  for  the  e  modes.  Engelking  et  al. 

[10]  and  lnoue  et  al.  [11]  have  published  experi¬ 
mental  values  for  several  frequencies.  By  combining 
theoretical  predictions  and  experimental  data  with 
observations  on  similar  molecules,  such  as  CH3F*,  we 
assign  the  following  values  to  the  vibrational  eigen¬ 
values  for  the  methoxy  radical  (directly  observed 
values  are  underlined:  3006  (e),  2930  (a,),  1380  (e), 
1325  (a,),  1060  (a,),  690  cm-' (e). 

These  frequencies  correspond  to  a  methoxy  radical 
zero-point  vibrational  energy  equal  to  22.1  kcal  mole-1. 

3.  Dissociation  energy  predictions 

Because  linked-diagram  derived  methods  are  size 
extensive  [2] ,  the  electronic  energy  change  for  a  dis¬ 


sociation  reaction,  A  -*■  B  +  C,  may  be  written 

AEel  =  Ee0)  +  Ecl(C)-Ea(A). 

Since  each  of  the  reactions  considered  here  contains 
the  hydrogen  atom  as  one  of  the  products,  AEei  re¬ 
quires  a  correlated  calculation  for  only  the  reactant 
and  the  molecular  product.  Table  2  contains  predicted 
electronic  energies  for  methanol,  CH3OH,  methoxy 
radical,  CH3O,  formaldehyde,  CH20,  formyl  radical, 
HCO,  carbon  monoxide,  CO,  and  the  SCF  energy  for 
the  hydrogen  atom.  Table  2  also  contains  the  zero- 
point  vibrational  energy  for  each  molecule  in  the 
study  [10-15). 

The  information  summarized  in  table  2  is  sufficient 
to  determine  the  zero-degree  heat  of  reaction  for 
each  of  the  four  reactions  included  in  this  study.  We 
summarize  the  values  for  each  of  the  molecules  in  the 
third  row  of  table  3.  For  most  molecules,  there  is  no 
experimental  data  available  to  compare  with  the  pre¬ 
dicted  AffR(0  K)  values.  An  exception  is  formaldehyde. 
Reilly  et  al.  [  16]  have  confirmed  that  the  threshold 
energy  for  the  dissociation  of  formaldehyde  to  radical 
products  is  86.0  ±  1.0  kcal  mole-1.  The  SDQ  MBPT(4) 
calculations  predict  that  the  threshold  energy  equals 
85.0  kcal  mole-1.  In  this  case,  the  agreement  between 
theory  and  experiment  is  very  good. 

In  addition  to  the  experimental  data  of  Reilly  et 
al.  [16],  experimental  values  for  room-temperature 
heats  of  reaction  are  available  for  methanol  [17] , 
formaldehyde  [18]  and  the  formyl  radical  [19],  To 
compare  the  theoretical  predictions  with  these  heat 
of  reaction  data,  we  must  estimate  the  effect  of  the 
rotational  and  translational  degrees  of  freedom,  and 
the  effect  of  thepK  work  term  must  be  included. 


Tabic  2 

Electronic  energy  predictions  and  zero-point  energies  for  CH,0  molecules 

Molecule  Electronic  energy  (hartrec)  Zero-point  energy 

-  (kcal  mole  1 ) 


SCF  (UHF) 

D  MBPT(4) 

SDQ  MBPT(4) 

CH30H 

-115.0744 

-115.4942 

-115.4920 

31.1  ( 12] 

CH30 

-114.4532 

-114.8248 

-114.8232 

22.1  [9] 

ch2o 

-113.8974 

-114.2832 

-114.2829 

16.2  (13) 

CHO 

-113.2769 

-113.6343 

-113.6363 

7.8  (14] 

CO 

-112.7665 

-113.1105 

-113.1124 

3.1  (15] 

H 

-0.49778 
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Table  3 

Dissociation  energies  predicted  by  SDQ  MBPT(4)  results 
(in  kcal  mole'1) 


CII3OH 

CH30 

CHjO 

CllO 

Ai*1 

107.3 

26.6 

93.4 

16.4 

Afvib 

-9.0 

-5.9 

-8.4 

-4.7 

A£r(0  K) 
Af.-U(300  K) 

98.3 

20.7 

85.0 

11.7 

0.9 

0.9 

0.9 

0.9 

Afrot(300  K) 

0 

0 

0 

-0.3 

A%(300  K) 

99.2 

21.6 

85.9 

12.3 

^R.300 

99.8 

22.2 

86.5 

12.9 

Translational  and  rotational  energy  contributions  are 
estimated  by  assuming  that  each  degree  of  freedom 
contributes \kT.  Thus,  for  T  =  300  K,  each  degree 
of  freedom  contributes  «=  0.3  kcal  mole-1.  The 
translational  and  rotational  energy  changes  for  each 
reaction  are  listed  in  table  3.  Combining  these  con¬ 
tributions  with  the  zero-degree  dissociation  energy, 
we  obtain  the  energy  of  the  reaction  for  T  =  300  K. 
Finally,  we  assume  that  pV=  An(RT),  and  add  this 
contribution  to  A/rR(300  K)  to  obtain  predicted 
values  for  the  heat  of  reaction,  AWR(300  K),  which 
are  listed  in  the  last  line  of  table  3. 

Table  4  contains  a  comparison  of  the  heats  of 
reaction  predicted  by  our  SDQ  MBPT(4)  calculations 
with  experimental  values.  Several  conclusions  may  be 
gleaned  from  the  theoretical  results  summarized  in 
table  4.  First,  the  inclusion  of  correlation  effects  is 
necessary  to  predict  dissociation  energies  with  chem¬ 
ical  accuracy.  Ab  initio  self-consistent-field  calcula¬ 
tions  predict  dissociation  energies  substantially  dif¬ 


ferent  from  experimental  values,  bach  of  the  corre¬ 
lated  treatments,  D  MBPT(4)  and  SDQ  MBPT14), 
predicts  dissociation  energies  that  agree  well  with 
experiment.  There  is  little  reason  to  prefer  the  more 
expensive  SDQ  MBPT(4)  calculations  to  the 
D  MBPT(4)  calculations  for  these  examples.  In  each 
case  reported  here,  the  contributions  of  the  single- 
and  quadruple-excitation  diagrams  at  fourth  order 
almost  cancel.  This  is  not  always  the  case,  so  there 
are  instances  where  SIX?  MBPT(4)  is  preferred.  The 
fourth-order  triple  excitations  can  also  be  important 
for  some  examples  [20] .  Finally,  in  each  case  where 
theoretical  predictions  can  be  compared  with  experi¬ 
mental  results,  the  predicted  dissociation  energy  is 
less  than  the  experimental  value.  This  is  not  sur¬ 
prising.  The  finite  basis  sets  used  in  molecular  cal¬ 
culations  describe  the  product  molecules  better  than 
they  describe  the  reactant.  This  characteristic  implies 
that  the  electronic  energy  difference  between  prod¬ 
ucts  and  reactant  usually  will  be  too  low.  Reliable 
estimates  of  the  magnitude  of  this  error  are  nui 
available,  but  our  experience  [21-24]  suggests  that 
dissociation  energies  predicted  by  MBIT,  for  breaking 
a  single  bond  in  a  polyatomic  molecule,  and  for  the 
quality  of  basis  set  used  in  these  calculations,  typi¬ 
cally  underestimate  experiment  by  1-3  kcal  mole-1. 
On  the  other  hand,  atomization  energies  are  much 
poorer  due  to  the  limited  basis  sets  [21]. 

4.  Heat  of  formation  of  the  methoxy  radical 

Two  experimental  values  for  the  heat  of  formation 


Table  4 

Comparison  of  predicted  and  experimental  dissociation  energies  (in  kcal  mole  1 ) 


Molecule 

Predicted  jqq 

Experiment 

SCF  (UHF) 

D  MBPT(4) 

SDQ  MBPT(4) 

tfR.300 

ref. 

CH3OH 

69.9 

100.2 

99.8 

100.9 

[10] 

104.0 

1171 

CII3O 

32.9 

23.1 

22.2 

cu2o 

70.1 

86.4 

86.5 

88.2  ±  1.6 

1181 

87.5  ±  1.0  a) 

[16] 

CllO 

4.4 

12.9 

12.9 

15.5  ±  1.5 

118] 

17.0  ±  2.0 

119] 

a)  Wc  estimate  this  value  based  upon  Ihe  jqq  value  of  86.0  ±  1  [16]. 
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CH3O 

CHO  + 

h2 

electronic  energy  al 

-114.8232 

-113.6363 

-1.1656 

AE*1 

(0.0213  a>) 

13.4  t» 

zero-point  energy  1,1 

22.1 

7.8 

6.2  [15] 

ae^p 

-8.1  b> 

AEr(OK) 

5.3  b> 

rotational  energy  (300  K) 

0.9 

0.9 

0.6 

ARrot(300  K) 

0.6  b> 

translational  energy  (300  K) 

0.9 

0.9 

0.9 

A£'r 

0.9  b> 

A£\(300  K) 

6.8  b> 

_  *4,300 

7.4  b> 

al  In  hartrec.  bl  In  kcal  mole  '. 


of  the  methoxy  radical  have  been  reported  recently  and  H2,  we  estimate  a  zero-point  energy  change  of 

[10,17],  Batt  and  McCulloch  [17]  derived  a  heat  of  —8.1  kcal  mole-1  for  the  reaction.  When  rotational 

formation  equal  to  3.9  kcal  mole-1  using  parameters  and  translational  energy  effects  and  our  estimate  of 

obtained  from  chemical  kinetic  data.  Engelking  et  al.  the  pV  terms  are  included,  we  obtain  a  value  for  the 

[10]  measured  the  electron  affinity  of  the  methoxy  heat  of  reaction  equal  to  7.4  kcal  mole-1.  Assuming 

radical,  and  estimated  the  heat  of  formation  for  that  AHf^oo  f°r  ,lle  f°rmyl  radical  is  9.0  kcal  mole-1 

methoxy  to  be  0.7  kcal  mole-1.  Since  our  calculations  [25] ,  we  predict  that  the  heat  of  formation  for  the 

predict  heats  of  reaction  for  two  processes  that  in-  methoxy  radical  is  1.6  kcal  mole-1, 

elude  the  methoxy  radical,  by  combining  these  re-  The  lack  of  agreement  among  the  predicted  values 

suits  with  published  heat  of  formation  data  we  can  for  the  heat  of  formation  is  not  surprising.  It  is  ap- 

obtain  a  heat  of  formation  for  CH3O.  We  assume  the  parent  from  the  results  of  this  work  and  the  results  of 

following  heats  of  formation:  CH3OH,  -48.0  kcal  theoretical  studies  using  alternative  methods  [1,3] 

mole-1;CH20,  -26.0  kcal  mole-1,  H,  52.1  kcal  that  theoretical  calculations  underestimate  dissocia- 

mole- 1  [25] .  Combining  the  dissociation  energy  of  tion  energies.  For  the  H-dissociation  reactions  re¬ 
methanol  with  the  heats  of  formation  of  methanol  ported  here,  and  in  other  work  done  using  these 

and  hydrogen,  we  predict  the  methoxy  heat  of  for-  methods  [21 ,26] ,  our  dissociation  energy  errors 

mation  to  be  -0.3  kcal  mole-1.  On  the  other  hand,  range  from  I  to  3  kcal  mole-1.  Considering  the  sign 

combining  the  methoxy  radical  dissociation  energy  of  this  error,  it  is  reasonable  to  suggest  that  the  heat 

with  heats  of  formation  for  formaldehyde  and  hy-  of  formation  of  the  inethoxy  radical  should  be  greater 

drogen,  we  predict  the  methoxy  heat  of  formation  than  0.3  kcal  mole-1  and  less  than  3.9  kcal  mole-1, 

to  be  3.9  kcal  mole-1.  The  heat  of  formation  obtained  by  studying  the  hy- 

ln  addition  to  the  two  hydrogen  elimination  reac-  drogen  molecule  elimination  reaction  of  CH3O  sug- 

tions  already  considered,  the  heat  of  formation  for  gests  a  heat  of  formation  for  the  radical  of  1 .6  kcal 

CH3O  can  be  obtained  by  considering  the  elimination  mole-1  that  lies  between  the  two  extreme  values, 
of  the  hydrogen  molecule  from  the  methoxy  radical.  Taking  the  likely  errors  into  account  we  suggest  a 

giving  the  formyl  radical  and  H2  as  products.  The  value  equal  to  2.0  ±  3.0  kcal  mole-1, 

pertinent  information  needed  to  estimate  the  heat  of 
reaction  for  this  process  is  summarized  in  table  5.  The 
results  of  SDQ  MBPT(4)  calculations  predict  a  classi-  References 
cal  dissociation  energy  Dc  equal  to  13.4  kcal  mole-1. 

Combining  the  estimated  frequencies  for  the  methoxy  1 1 1  J.A.  Pople,  R.  Krishnan,  H.B.  Schlcgel  and  J.S.  Binkley, 

radical  with  empirical  values  for  the  formyl  radical  Inicrn.  J.  Quantum  Chem.  14  (1978)  545. 
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